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ABSTRACT 



Aims. We study the physics of cyclotron line formation in the high-energy spectra of accreting X-ray pulsars. In particular, we link 
numerical predictions for the line profiles to results from observational data analysis. Therefore, first we investigate the theoretical 
predictions and the significance of our model parameters, and second we aim at the development of a model to fit cyclotron lines in 
I . observational data. 

Methods. Simulations were performed using Monte Carlo methods. The data were extracted with HEADAS 6.1.1 and INTEGRAL 
OSA 5.1. A convolution model for the cyclotron line shapes was implemented for the XSPEC spectral analysis software package and 
CZ3 , for data packages compatible with XSPEC local models. 

Results. We predict the shapes of cyclotron lines for different prescribed physical settings. The calculations assume that the line- 
forming region is a low-density electron plasma, which is of cylindrical or slab geometry and which is exposed to a uniform, sub- 
critical magnetic field. We investigate the dependence of the shape of the fundamental line on angle, geometry, optical depth and 
temperature. We also discuss variations of the line ratios for non-uniform magnetic fields. We have developed a new convolution and 
interpolation model to simulate line features regardless of any a priori assumed shape of the neutron star continuum. Fitting RXTE and 
INTEGRAL data of the accreting X-ray pulsar V0332+53 with this model gives a qualitative description of the data. Strong emission 
wings of the fundamental cyclotron feature as predicted by internally irradiated plasma geometries are in principle observable by 
todays instruments but are not formed in V0332+53, hinting at a bottom illuminated slab geometry for line formation. 
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■ 1. Introduction served cyclotron line shapes are energetically resolved in de- 

>_ , _ . ^ „ , tail for several sources, thus demanding an in-depth study of 

Cyclotron Resonance Scattering Features (hereafter called • c .■ A , .. , . c • , ,. 

• i— i . _£,„_, . , „ , ,. „. n ,. ,. , their formation. A better understanding of their complex line 

, CRSFs or simply cyclotron lines"), first discovered jn toe s pec- sh wouM reyeal much about ^ h ^ ^ whkh 

■ trumofHerX-1 dTrumper et aU | l977 fc| Trumper et al| | l978|) , are driyes ^ fascinati but t0 tod onl , understood 
Ctf ; observed as absorption lin es in the spec traofhighly magnetized ses of accretion in x lsars At nt> the ob _ 

accreting X-ray pulsars jHeindl et al. [ £004j). Their scientific seryed ^ ^ usually modeled fe Gaussian Qr Lorentzian 
importance lies in providing the only direct method currently gh (Mi hara1iiLl[l990; MakishimalULl \\9M). For mod- 
known f or the determina t^f the magnetic field ofaneutro" e ling the fundamental feature, one must sometimes use several 
star (e.g. | Harding & Lai EQQ6J; l Orlandini ^& Fiumei |200l|) The Lorentzians or Gaussians to obtain an acceptable fit of its non- 
me energy E cyc of the fundamental CRSF is approximately re- sh ( Kj-evkenboh^iliLl [2005). 
lated to the magnetic field strength of the star by the 12-B-12 

rule" In this paper, using Monte C arlo simulations based o n an 

E cyc ~ 11.57 keV • An, (1) improved version of the code of lArava & Harding! (H~999) and 

lArava-Gochez & H arding (2000), we calculate line features for 

where B n denotes the m agnetic field in units of 10 12 Gauss typical neutron star spectra and infer the line profiles under the 

dCanuto & Venturat |1977|). Due to the gravitational redshift, assumption of physical parameters such as the accretion geome- 

Z, the observed line energy is shifted by a factor of try, the viewing angle, and the plasma temperature. Based on our 

1/(1 + z) with respect to £ cyc . Today, with the access to simulations, we develop a new local XSPEC model, cyclomc, 

data from satellites as BeppoSAX, the Rossi X-ray Timing for cyclotron lines and show first results from fitting observa- 

Explorer (RXTE), the International Astrophysics Gamma-Ray tional data with cyclomc. 

Observatory (INTEGRAL), and Suzaku, the diagnostic poten- ,_. 

tial of CRSFs has grown anew: with these instruments the ob- The outllne of the P resent work 18 as follows: In Sect. |2] we 
summarize the basic theory of cyclotron line formation and de- 
Send offprint requests to: G. Schonherr, scribe standard scenarios for the accretion process onto the neu- 
e-mail: gschoen@astro.uni-tuebingen.de tron star. We report key results from past observations and give 
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an overview of different numerical approaches to modeling cy- 
clotron lines for neutron stars. Sect. [3] focuses on the modeling 
approach taken here. We try to motivate our steps from the point 
of view of an observer, however, some rather technical discus- 
sion is necessarily included in that section. Theoretical predic- 
tions from Sect. |4]are complemented in Sect.[5]by a comparison 
of the model with real data. Finally, in Sect. [6] we summarize our 
results and discuss future steps for cyclotron line modeling. 

2. Overview 

2. 1 . CRSF formation in X-ray pulsar spectra 

Cyclotron lines are found in the spectra of accreting neutron 
stars in binary systems with magne tic field strengths of t he or- 
der B ~ 10 n -10 13 Gauss. See, e.g jHarding & Lail d2006l) for a 
review of the physics of strongly magnetized neutron stars. X- 
ray pulsars are thought to have masses of M ~ 1 .4 Mq, radii 
R ~ 10 6 cm, and luminosities of 10 34 — 10 38 ergs -1 . They accrete 
matter from a strong stellar wind or by Roche lobe ove rflow, 
usually through an accretion disk (Ghosh & Lambl |l978i) . The 
magnetic field of the neutron star disrupts the flow of matter at 
the Alfven radius, where the magnetic field pressure equals the 
ram pressure of the flow, and the mat ter is funneled along th e 
field lines onto the magnetic poles (Basko & Sunvaev, 1976), 
reaching free-fall velocities of ~ 0.4 c. When discussing pro- 
cesses of accretion in the following, we restrict ourselves to the 
field-dominated volume near the neutron star surface at the mag- 
netic poles. Inverse Comptonization of soft photons in the de- 
celerated plasma produces photons in the X- and gamma-ray 
regime. The emission characteristics of this radiation depend on 
the mass accretion rate, M. For large M, a shock front develops 
at some distance from the neutron star surface, which does not 
permit the upscattered photons to escape vertically from the ac- 
cretion column, i.e., parallel to the B-field. As a result, a 'fan 
beam' emission pattern forms (Fig. [TJ. As was first shown by 
Bask o & Sunvaevl dl976l) . the critical luminosity for shock for- 
mation, L*, is 



L* = 2.72 • 10 



37 



CTj 



V°"ll°"-L 



(r \( M 
\r)\Mq 



erg s 



(2) 



where ro is the polar cap radius, <x T is the Thomson scattering 
cross section and cry and cr ± are the energy averaged cross sec- 
tions for the scattering of photons which propa gate in paralle l 
and perpendicular to the magnetic field direction (Becker, 1998). 
For small M, on the other hand, i.e., for L < L* , the radiation 
is emitted from an accretion mound such that most photons are 
emitted parallel to the B-field in a 'pencil beam' pattern. 

Before emerging from the line-forming region, the high- 
energy photons undergo scattering processes with the electrons 
in the relativistic plasma of the accretion column. The scatter- 
ing cross section is resonant at energies determined by the sep- 
aration of the Landau energies, the discrete energy levels of the 
electrons: when the strength of the magnetic field B approaches 
the critical field strength, B clit = (m 2 c 3 )/(efi) = 44.14 • 10 12 G, 
the de Broglie radius of a plasma electron becomes compa- 
rable to its Larmor radius. Quantum mechanical treatment of 
the electrons' motion perpendicular to t he magnetic field lines 
dMeszarosL [l99l iDau ghertv & Harding, QHll) reveals a quanti- 
zation of the electrons' momenta p±/(m e c) = n{B/B cl i t ). This 
translates into discrete energy levels, where the fundamental 
Landau level is given by the Yl-B-Yl rule (03 and the higher 
harmonics have n times this energy. For photon-electron scatter- 
ing, relativistic effects lead to a slightly anharmonic spacing of 



the resonant photon energies. Due to the large scattering cross 
section at the resonances and due to the quasi-harmonic spac- 
ing of the thermally broadened Landau levels, photons of ener- 
gies close to the Landau level energies may not escape the line- 
forming region unless inelastic scattering has slightly changed 
their energy. Consequently, absorption features in the photon 
spectrum are observed at 



E„ — m e c' 



1 + 2n(B/B crit ) sin 2 0-1 1 



sin 2 6» 



1 



(3) 



where m e is the electron rest mass, c the speed of light, 6 the 
angle between the incident photon direction and the magnetic 
field vector, and z is the gravitational redshift at the radius of the 
line-forming region. Note that we enumerate the cyclotron lines 
starting at n — 1 , and refer to them as the first or fundamental line 
at the energy E^ c = E\, followed by the second, third, fourth, 
etc. harmonics (n — 2, 3, 4, • ■ ■ ). The gravitational redshift at the 
neutron star surface is approximate!^ 



1 



1 



(4) 



1GM 
Rc 2 



which gives z ~ 0.3 assuming the typical neutron star param- 
eters given above. The thermal motion of the electrons parallel 
to the magnetic field lines remains free and is characterized by 
the parallel electron temperature, T e , which will be introduced 
in Sect. 133] 

2.2. Observations 

In 1976, the first cyclotron line was det ected in the X-ray spec- 
trum of Her X-l (Tri imper et all 119771) . First interpreted as an 
emission feature at 53 keV, the line was lat er proposed to be in 
absorption with theoretical arguments by iNagell d 198 lb . Since 
the discovery of the Her X-l cyclotron lin e, more sour c es ex- 
hibiting CRSF s have been observed (e.g., iHeindl et al. 



2004; 



IStaubertll2003HCoburn etail. 120021: ISantangelo et al.ll2000t) . At 
the time of writing, more than 16 accreting pulsars with securely 
detected cyclotron lines with surface magnetic fields in the range 
of 1-5 • 10 12 G were known. The record holder with respect 
to the number of lines detec ted is 4U 0115+63 dHeindl et all 
1 19991: ISantangelo et all fl 999) where five harmonics were found 
(IHeindl et al.Ll2000l) . 

CRSF sources are regular targets for observations. The 
progress over the last decades in energy resolution of instru- 
ments on satellites like BeppoSAX, RXTE, INTEGRAL, and 
Suzaku has led to excellent observational data of many interest- 
ing objects with complex cyclotron line features and has made 
high-quality phase resolved spectroscopy possible. As a result, 
many interesting characteristics of CRSFs are known today, 
awaiting a deeper explanation than given by the simple picture of 
line formation which was outlined in the previous section. Some 
key results from observational studies are: 

1. The profile of the fundamental line is resolved, is clearly 
non-Gaussian, and exhibits a complex shape. 

2. The second harmonic generally appears deeper than the fun- 
damental line. 

3. Significant variations of the line parameters of the CRSFs 
with the pulse phase are observed for some sources. 



1 Strictly spoken, Eq. Q is exact only for a non-rotating, spherically 
symmetric, uncharged mass. 
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4. The line ratios are not necessarily harmonic. The deviations 
from the harmonic energies in some spectra are too large to 
be explained only by the basic relativistic corrections im- 
plied by Eq. ||3). 

5. The line position of the funda mental CRSF can vary with the 
sourc e luminosity. Negative jMiharal 1995 : Mowlavi et al.L 
|2006j; Eakaiimaet al., 2006; Tsvgan kov et al l 120061) and 
positive dLa Barbera et al.Ll2005LIStaubert et alll2007l) linear 
energy-to-luminosity correlations have been found. 

2.3. Numerical models 

There are two very different approaches to modeling the radia- 
tive transfer in the accretion column: solving finite difference 
equations and Monte Carlo simulations. 

2.3.1. Solving difference equations 



Motivate d by the Her X-l l i ne det ection, iNagell (Il980l Il98ll) 
and later Meszaros & NageJ d 1985b employed Feautrier meth- 
ods in order to solve the radiation transfer equation. They per- 
formed two sets of calculations, treating effects of anisotropy 
and Comptonization separately. Having first presented a line- 
forma tion mechanism for a cyclotron emission feature (Nagel, 
1980), in a later paper considering Comptonization effects, 



Nagel (119811) then favored the Her X-l line to appear in absorp- 



tion. For the combi n ed ef fects of anisotropy and Comptonization 
Mesz aros & Na gel (1985) compared model predictions for dif- 
ferent geometries (slab and cylinder geometry with internal or 
external illumination) and discussed variations with the angle of 
the emergent spectra. Their ap proach was later refined by the in- 
clusion of higher harmonics ([Alexander & MeszarosL fl99lh or 
by including radiati on pressure and temp erature corrections in 
the atmosphere dBulik et al.L 1 19921 1 1995b . Recently, the influ- 
ence of a non-uniform magnetic field in the line-forming region 
on the formation of CRSFs has been investigated with similar 
techniques by Nishimura, who found a varia tion of the line ra- 
tios of the CRSFs for the case of a dipolar (iNishimural 120031) 
and for a linearly varying (Nishimura, 2005) magnetic field. 



where the slab normal ma y have any direction wi th respect to the 
B-field were discussed by Isenberg et al. ( 1998a b). For the case 
where the slab normal is perpendicular to the magnetic field vec- 
tor, their results are comparab le to assuming a cyl inder geometry 
for the line-forming region. Isenb erg" et alj dl998bl) distinguished 
between line shapes of optically thin and optically thick matter. 
As one key result these authors found that the line wings dis- 
appear either for the 1-1 geometry and optically thick media or 
for the 1-0 geometry and optically thin media. However, none 
of these scenarios could explain the observed fundamental shal- 
low and broad features due to the high equivalent width of the 
fundamental in both cases. 

Inspired by the detection of up t o two cyclotron linefl 
at 50 and lOOkeV dKendziorra et all 11994 and at HOkeV 
dGrove et all 1 19951) during subsequent outbursts of the transient 
source A0535+26 in 1989 and 1994, Araya and Harding pre- 
sented a new set of Monte Carlo simulations for very hard spec- 
tra of X-ray pulsars with near-critical fields dArava & Harding! 



ray 1 

1996. ll999tlArava-G6chez & Hard ing. 200$. For a low-density 
plasma and hence low continuum optical depths, they produced 
spectra for slab ( 1 - 1 ) and cylinder geometry of a plasma threaded 
by near-critical magnetic fields and discussed the influence of 
parameters as geometry, optical depth and anisotropy of the pho- 
ton source on the line shapes. The results presented in this paper 
are based on their approach. 



3. Modeling CRSFs 

3.1. Aims 

The key objective of this work is to obtain a physically motivated 
model for CRSF formation which is directly comparable to ob- 
servational data. Firstly, such a comparison is fundamental when 
testing and reconsidering the validity of the model. Secondly, 
once the model has reached a well-developed state, only its sim- 
ple applicability to real observational data provides the means 
for a systematic investigation of CRSF sources. 

3.2. Methods 



2.3.2. Monte Carlo simulations 

lYahel (fl979b was the first to use Monte Carlo simulations for 
simulating the CRSF formation in the atmosphere of a mag- 
netized neutron star. He considered the formation of pulse 
profiles and X-ray spectra and found that the Her X-l fea- 
ture could indeed be reproduced as a consequence of reso- 
nant scatterings of extraordin ary polarized photons. Two years 
later, Pravdo & Bussardl d!981l) calculated angle-dependent pul- 
sar spectra, including relativistic corrections to the Compton 
cross section and considering polarization dependence. Focusing 
on the continuum spectral shape they f ound a hard e ning o f the 
spectra towards the magnetic equator. IWang et ail (Q989) per- 
formed Monte Carlo simulations for the geometry of a plane 
parallel slab with the slab normal parallel to the B-field vec- 
tor, and the plasma being ill uminated from below . This geometry 
was na med '1-0 geometry' dFreeman et allll99"9l : llsenberg et al.L 
1998b]), in contrast to the '1-1 geometry' which indicates a ho- 
mogeneous slab illuminated at the mid plane. While in the 1-0 
geometry photons which return to the source plane after scatter- 
ing are absorbed ('reflected photon flux'), in the 1-1 geometry 
photons may cross the source plane and the reflected and trans- 
mitted flux are symmetric. Results from a generalized model, 



In order to achieve the desired flexibility, we base our model 
on Monte Carlo simulations using a revised, generalized version 
of the Araya & Harding code. A new key feature is a Green's 
functions approach giving independence from any continuum 
model assumed. Araya & Harding's implementation of an in- 
ternally irradiated slab geometry is generalized to include also 
the case of illumination from the bottom (1-0 geometry). More 
information on the Monte Carlo implementation and details of 
the Green's function approach are given in Sect. 13.41 Due to the 
variety of sources and the uncertainty of the general physical 
picture, calculations are performed on a large multidimensional 
parameter grid. All simulation results are merged into archives 
in the form of FITS tables which are available from the authors. 
Line features for X-ray pulsar spectra for different physical set- 
tings (as outlined in the next section) within our parameter scope 
may be produced from these tables with a special convolution 
and interp olation model, also implemented as a local model for 
XSP EC dArnaudL 1 19961) and o ther analysis packages such as 
ISIS dHouck & DenicolaLl2000b . 



2 Re cently it was shown that both lines cla imed by[ Kendziorra et ap 
d!994f) are present in the source spectrum dKretschmar et all 120051 
IWilson & Fingel l2005: Ino ue et all 12001 ; ICaballero et aUl20d% . 
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3.3. Physical setting 

We simulate the propagation of photons through a medium of 
prescribed physical conditions. The photons interact with the 
electrons in the medium via resonant scattering processes. The 
conditions in the line-forming region of the electron plasma ar e 
governed by the following parameters (Araya & Harding, 1999): 

1 . Magnetic field [B] As a first approximation we consider a 
neutron star with a magnetic field which is assumed to be 
uniform on the scale of the line-forming region. The field 
strengths simulated are between 1 • 10 12 and 7 • 10 12 Gauss, 
encompassing the whole range of fi-fields found in cyclotron 
line sources. In Sect. 14.31 a possible generalization to non- 
uniform magnetic fields is described. 

2. Plasma electrons [T e , f(p e ), n] We consider a low-density 
thermal plasma. We also assume that all electrons are ini- 
tially in their fundamental Landau state n — 0. This assump- 
tion is justified by the very high cyclotron radiative decay 
rate for sub-critical fields 



r^S-lO 15 ^" 1 
compared to the collisional excitation rate 

r col = 5 ■ 10 8 (» e /10 21 cm- 3 )B:, 3/2 s- 1 



(5) 



(6) 



dLatall 119861: iBonazzola et all Il979h . For their motion par- 
allel to the Z?-field vector, we assume a thermal distribution 
of the electrons with their parallel momenta p e given by a 
relativistic Maxwellian distribution 



f(p e )dp e oc exp 



m e (r 



2 \\ 
1 



kT e 



dp e 



(7) 



where T e is the parallel electron temperature and k is the 
Boltzmann constant. In the literature, T e is often linked 
to th e strength of the ma gnetic field (e.g. Lamb et al. . 
1 9901: llsenberg et al.L 1 1 998bHArava & Hardingl 1 1 996L 1 1 999 : 
Araya-Gochez & Harding[|2000l) . We take these studies into 
account in order to determine the order of magnitude of the 
plasma temperature, but leave T e as a free parameter in our 
simulations in order to keep the model's flexibility (see also 
Sect. 14.41 ). We also assume a slowly sinking plasma where 
bulk plasma motion may be neglected. 
Optical depth [tj] The Thomson optical depth tj of the 
plasma is prescribed. The optical depth for cyclotron scatter- 
ing r cyc relates to the Thomson optical depth as 



0~cyc 
^"cyc — > 



(8) 



implying a scattering optical depth T cyc which can be a fac- 
tor of ~10 5 larger than tt at the resonances of cr cyc . Fig. [2] 
shows the thermally averaged cross section (cr^/crj as a 
function of e nergy and an gle, calculated as a second order 
QED process (ISinal [19961) . Besides a highly resonant behav- 
ior of the cross section at the Landau energies, Fig. [2] also 
illustrates the angle-dependent relativistic shift in the reso- 
nances as well as the thermal broadening of the profiles. We 
calculate line features for Thomso n optical d epths be twee n 
Tt — 1 ■ 1(T 4 and t t = 3 • 1(T 3 dArava & Harding! [1999). 
Depending on the plasma geometry, the resulting mean free 
path of a photon in the line-forming region is different for 



the same trajectory. The Thomson optical depths we simu- 
late correspond to electron column densities N e = Tt/o"t 
between 1.5 ■ 10 20 cirT 2 and 4.5 • 10 21 cirr 2 . These values 
are comparable to values assumed in other recent numerical 
or analytical studies. For instan ce, column densities of jV e ~ 
10 21 -10 22 crrT 2 are assumed bv lNishimura l d2003ll2005l). and 
values of ~ 10 22 crrT 2 are inferred by iBecker & Wolfll 
( 120071) for the sources Her X-l, LMC X-4, and Cen X-3. 
Geometry [geo: sl/cy] We distinguish two basic geome- 
tries of the line-forming region motivated by the comple- 
mentary 'standard' pictur es of accretion depicted in Fig. Q] 
(see also section 12.11 and [Basko & Sunvaevt Il976l) : for the 
case of flow stopping through nuclear collisions at the 
surfa ce, we adopt the geometry of a thin, pla ne-parallel 
slab (Mes zaros et all 119831: iHarding et all fl984). Radiative 
shocks or shocks from collisionless instabilities, on the other 
hand, require a cylindr i cal sha pe of the X-ray emitting region 
(e.g., IBecker & Woil 120071 and references therein). The 
heights of slab atmo s phere s are expected to b e significantly 
smaller (lLamb et all 1 1973b IWang et all 1 19881) than typical 
radii, h <g rp, of the accretion mound dOstriker & Dav idson. 
19731: IBecker & Wolfll [2007). For cylindrical geometries, 



Becker & Wolff (2007) investigated a more complicated ra- 



diative shock structure with a velocity gradient. For several 
X-ray pulsars, these authors found height to width ratios, 
Zmax/ro, of the emitting region of the order 10 2 , confined 
to a volume smaller than the one of trapped radiation. Our 
simulations also suggest that practically all photons escape 
before reaching such boundaries, justifying our assumption 
of infinite extended plane-parallel line-forming regions. For 
the cylinder, we assume that the plasma is irradiated inter- 
nally with the photon source located at the cylinder axi s 
dArava & Hardingl [l999; Araya -Gochez & Harding, 2000). 
Two locations of the source plane are considered for the slab. 
The scenario of a line formation region above an isotropi- 
cally emitting source is realized by a bo t tom-illuminated slab 
(Free man et all 119991: llsenberg et all Il998bt IWang et all 
I1989I) . while placing the photon source at the midplane 
of a plane-parallel slab is representative of the scenario of 
line formation in an isothermal, semi-infinite atmosphere 
dSlater et al.L[l982h . 

The optical depth the photons see into direction Q when cov- 
ering an optical depth Atj along the slab normal or perpen- 
dicular to the cylinder depends on the geometry as: 



Ar T (6») = 
At t (6») = 



At t 

sin# 
Ar T 

cos 6 



(cylinder geometry) (9) 
(slab geometry) . (10) 



3.4. Technical realization 



We use a Monte Carlo method to simulate the reso- 
nant scattering processes between incident photons and 
plasma electrons which lead to the formation of cyclotron 
line features. The code use d is a revised v e rsion of the 
code of Araya & Harding dArava & Hardingl Il996l 119991: 
lArava-Gochez & Hardingl |2000|) . The calculation of the rela- 
tivistic cross sections is done with a separate code by ISinal 
d!996l) . Resonant scattering with electrons up to the fourth har- 
monic is included. 

Besides relaxed geometrical constraints on photon injection, 
and technical modifications such as an improved angular and en- 
ergy resolution and increased statistics, the main difference of 
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"fan beam" 



"pencil beam" 



free falling plasma 



shock 




layer of radiation 




Fig. 1. Accretion geometries and radiation patterns. Left: "fan 
beam"/cylinder geometry. Right: "pencil beam"/slab geometry. 



our program with respect to the preceeding one is the Green's 
functions approach: In each Monte Carlo run we insert 10000 
photons of the same incident energy, E- m , pick for each photon a 
random angle 9 m (cos(f5i„) e (- 1 , 1)) with respect to the magnetic 
field direction, propagate them through the plasma (see below), 
and calculate from all final states the probabilities for photon 
redistribution into different energy and angular bins. The initial 
angular distribution of the photons is assumed to be isotropic. 
The relevant energy range \E m } for cyclotron line formation is 
assessed as follows: using the 12-5-12 rule and assuming quasi- 
harmonic spacing of the cyclotron lines, the energy range con- 
taining the first four Landau levels can be fixed independently 
from the magnetic field strength in terms of E/B12. Due to the 
link of resonant energy and magnetic field, the choice of this 
scale is important. First of all, it gives an optimized resolution of 
the CRSF features in the same way for different magnetic fields. 
Even more important, the choice of this scale is fundamental for 
later interpolation of the Green's functions: as the resonant en- 
ergies are directly linked to the magnetic field, an interpolation 
of the line shapes in E/Bn-space ensures consistent results. For 
£in/£i2/keV G [6,48], we obtain a grid of Green's functions 
G(E m — > E out , 9 out ). E[ n is sampled by 161 Monte Carlo runs, 
the resolution of the redistributed energies E out is given by an 
internal energy binning of 640 bins. 

Each Monte Carlo photon is injected into the plasma with 
energy E m and direction 9 m . The photon is then propagated ac- 
cording to its mean free path l/(n e (o~{E m , cos #;„)) ) and an elec- 
tron is picked as a scattering partner. The electron is character- 
ized by its parallel momentum, p e (Eq.|7]i, and its Landau state n. 
According to the scattering cross section (obtained from interpo- 
lation of previously calculated and tabulated values as a function 
of B), the state of the electron-photon pair changes from its in- 
cident configuration (£' (0) , f3 <0) ) + {pf\n (0 ^) to a different state 
(E^ l \ 6> (1) ) + (pi , n (1> ). The new mean free path of the photon 
is calculated and the photon is propagated further. If the elec- 




20 40 60 80 100 

Energy [keV] 

Fig. 2. Cross sections for B/B cr [ t = 0.05 and kT e = 3keV, ther- 
mally averaged over the sampled electron momenta, and in units 
of the Thompson cross section cr T . The resulting profiles are 
shown for different angles 9 of the photon's direction with re- 
spect to the magnetic field vector. Solid lines: cos 9 = 0.005 and 
cost? = 0.875. Dotted, dashed-dotted and dashed lines (sharper 
towards broader): cosfl = 0.200,0.375,0.625. Figure based on 
lArava & Hardingl (Il999l) . 



tron remains in an excited Landau state n' > after scattering, 
another photon is emitted with (£ (2) , f9 (2) ) + (/?g 2) ,« (2) ) and pro- 
cessed further. This photon spawning can produce up to three 
secondary photons. Once a photon has escaped from the plasma, 
its contribution to the output spectrum is stored. Fixing the input 
angular distribution of the incident photons to be isotropic in the 
scope of the present work is done for reasons of simplicity, and 
in order to keep the computational expenses reasonable. In a fu- 
ture paper we will discuss the generalization to arbitrary angular 
photon distributions. Polarization of the photons is not included, 
however, polarized photons should yield a comparable picture 
(Wang Tt~a"Qll988l) for the low-density regime chosen. 

Calculations were performed in six-dimensional parameter 
space for a non-regular grid of points [B/B a -i t , T e ,p, tj, E m , geo] 
within the ranges B/B alt e [0.03, 0.15], kT e e [2.5, 20] keV, ji = 
cos(0 out ) e (0, 1), r T e [1 • 10 4 ,3-10 3 ], E m /B l2 /teV e [6,48]), 
and for slab 1-1, slab 1-0, and cylinder geometry respectively. 
The current parameter grid is resolved into (NB/B ah X Nt c xN^X 
N TT xN Ein xN seo ) = (16x4x8x4xl61x3)grid points requiring a 
simulation time of the order of 10 CPU hours on 2 GHz work- 
stations. The resolution was chosen such that the variation of 
the Green's functions between two points is sufficiently small to 
allow for interpolation. Hence, we can predict CRSFs by convo- 
lution of a continuum spectrum for any parameter combination 
(B/B cl -it, r e ,yu, tj, geo) on this grid as follows: First, we obtain 
the corresponding Green's functions G* by linear interpolation 
in all parameters except the geometry onE/B^. A set of Green's 
functions {G*(Ej — > Ej, 9k)} j for a fixed example input energy Ej 
and fixed physical setting is shown in Fig. [3] Second, we cal- 
culate the emergent photon flux F em (Ej, 9k), i.e., the number of 



6 



G. Schonherr et al.: A model for cyclotron resonance scattering features 




10 15 20 25 30 

E/B 12 [keV] 

Fig. 3. Green's functions plotted for one input energy 
Ej/B^/keV = 31.6 (dashed lines) against all output energies 
{Ej}j, and into four (of eight) cos 9 bins (cos 9 e [0.125, 0.250), 
[0.375,0.500), [0.625,0.750), [0.875, 1.0) from top to bottom). 
The physical setting is as follows: geo:sl, B/B CI [ t = 0.06, kT e = 
5 keV, tj = 3 • 10~ 3 . Ej corresponds to photon input at the third 
harmonic line. Most photons are redistributed by photon spawn- 
ing to the wings of the fundamental line. 



photons per ke V in the / h energy bin and k th angular bin (binned 
in cos 9), as a function of the (isotropic) incident continuum flux 
^cont convolving F cont with this interpolated set of Green's 
functions 

y, G*(£,- — > Ej, 9k) F cont (Ei)AEi 

F em (E h 8 k ) = — — J — — -. (11) 

AEj 

Note that because of this approach CRSFs for arbitrary contin- 
uum shapes can be calculated without rerunning the simulations. 



3.5. X-ray pulsar continua 

In the following, we show cyclotron resonance scattering fea- 
tures in folded full spectra. As continuum input, we chose an 
exponentially cutoff power law of the form 



F(E) = A-E- a exp(-E/E M d) 



(12) 



where A is the power law normalization, a the photon index 
and Efoid the folding energy. This spectral shape is the most 
simple phenomenological model which qualitatively describes 
X-ray pulsar spectra and is used here for illustrative purposes 
only. When considering real observational data, however, more 
complex con tinuum models like, e.g. a power law with a Fermi- 
Dirac cutoff dTanakaUl986l) or a ne gative and pos itive powerlaw 
with exponential cutoff ('NPEX') tiMiharal fl995l) can be used. 
Only recently, an analytical derivation of the spectral sh a pes of 
X-ray pulsar spectra was presented by iBecker & Wolff! (2005, 
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Fig. 4. Line profiles for cylinder geometry as a function of angle 
and optical depth. Left: continuum optical depth tj = 3 • 10~ 4 ; 
right: tj — 3 ■ 10~ 3 . The emergent spectra [photons ster -1 
s _1 keV -1 ] from all eight angular bins are shown (bottom to 
top: // = cosfl e [0.000,0.125), [0.125,0.250), [0.250,0.375), 
[0.375,0.500), [0.500,0.625), [0.625,0.750), [0.750,0.875), 
[0.875,0.1000)). In both panels, B/B ait = 0.05 and kT e = 
3.0 keV. The continuum photon flux is assumed to have a power 
law distribution with photon index a — 2.0 and an exponential 
high energy cutoff at the folding energy £f id = 40 keV. 



4. Results 

In this section we make theoretical predictions from our Monte 
Carlo simulations and discuss their implications on observed 
properties of cyclotron lines. Special emphasis is placed on the 
study of the line parameters, i.e., line position, line width and 
line depth, of the fundamental CRSF. 



4.1. Geometry and optical depth 



To c lose the link to the publicat i ons bylAraya & Harding! (Il999h 
and lArava-Gochez & Harding! d2000l) . we illustrate in Figs. |4] 
and [5] full spectra folded with the convolution model for a cho- 
sen physical setting in a similar fashion as shown in their work: 
spectra are depicted for slab and cylinder geometries, for two 
different values of tj and for different angular bins. The con- 
tinuum parameters are fixed to E{ n - 40keV and a — 2. The 
choice of the folding energy describes a rather hard continuum, 
which leads to a pronounced contribution from high-energy pho- 
ton spawning to the CRSFs (compare Sect. I4.5I >. Our results are 
in agre ement with the previous simulations of lArava & Ha rding 
( 1999). As an improvement with respect to the earlier results, we 
calculated the line shapes on a grid of 640 energy bins instead of 
80 and we resolve the depicted line shapes for eight angular bins 
of cos 9 instead of four. The improved statistics of effectively 
1.6 ■ 10 6 Monte Carlo photons per folded spectrum (increased 
from 5 ■ 10 4 ) gives us a well resolved picture of the line shapes, 
however, the centroid of the second harmonic would require even 
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Fig. 5. Same as Fig.|4]for slab 1-1 geometry. 
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Fig. 6. Same as Fig.|4]for slab 1-0 geometry. 



better statistics for continuum optical depths of tj ~ 3 • 10~ 3 . 
Although the second harmonic is not fully described in its core, 
the calculations clearly show that the second harmonic is more 
pronounced than the fundamental one, which has a more com- 
plex, broad and shallow shape, often with emission wings, which 
are strongest for the internally illuminated plasmas. The strength 
of those emission wings is expected as a consequence of the 
source photon location at the mid- plane of a slab and mid -axis 
of a cylinder. Fo r slab geometry, Isenberg et al. ( 1998b]) and 
iNishimural d2005l) showed comparisons of the slab 1-1 and the 
slab 1-0 geometry, and pointed out that these wings can be un- 



Fig. 7. Angular redistribution of the photons. For isotropic pho- 
ton injection, the percentage of the emitted photon flux into eight 
angular bins for the spectra in Figs. H] [5] and [6] is shown. Top: 
cylinder geometry, middle: slab 1-1 geometry, bottom: slab 1-0 
geometry. Solid lines: tj — 3 ■ 10~ 3 , dashed lines: tt = 3 ■ 10~ 4 . 



derstood as the remnants of a strong emission feature forming 
near the mid-plane. For very small optical depths, tj = 10~ 4 , we 
also find this feature in our simulations. Observations of sources 
with CRSFs have not been seen to exhibit strong emission wings. 
In Sect. 15.21 we will comment further on the observability of 
those features with modern instrumentation. 



4.2. Angular distribution 

Photons are injected isotropically into 20 angular bins in cos 8. 
Although the distribution of the initial photon directions is 
isotropic, a high degree of anisotropy arises after the pho- 
tons have been propagated through the plasma due to a highly 
anisotropic scattering cross section. Fig.|7]shows the angular re- 
distribution of the photons for all considered geometries and for 
different values of the optical depth. Internally the code keeps 
track of twenty cos 8 bins. We show the angular redistribution re- 
solved into eight final bins. For tj — 3 • 10~ 4 we observe a trend 
of an overall redistribution towards smaller 8 for slab geometry 
and a reverse trend for cylinder geometry. For a larger optical 
depth, tt = 3 • 10~ 3 , these trends increase for slab and decrease 
for cylinder geometry, where the curve in Fig.Qfiattens. This can 
be understood from the dependence of the scattering cross sec- 
tions on the angle (see Fig. |2j which implies that there is general 
trend of a photon redistribution by scattering towards smaller an- 
gles, i.e., larger cos 8, regardless of the geometry. The larger the 
optical depth a photon must pass, the more scatterings take place 
and the more dominant this effect becomes. This condition can 
also account for a less prominent trend in the slab 1-0 geome- 
try compared to the 1-1 geometry, as photons which experience 
many scatters by various crossings of the source plane are ther- 
malized, biasing the emerging radiation. For a fixed optical depth 
tj = 1 • 10~ 3 similar plots are shown resolved into all twenty bins 
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Fig. 8. Variation of the line parameters of the fundamental CRSF 
with the angle of emission. Solid lines: cylinder geometry, 
dashed lines: slab 1-1 geometry, dash-dotted lines: slab 1-0 ge- 
ometry. All values are obtained from fits of the line shapes de- 
picted in Figs. H] |5] and [6] and for a Thompson optical depth of 
r T = 3 • 1(T 3 . 



bv lArava-Gochez & Harding d2000t) . For r T ~ 8 ■ 1(T 4 , the an- 
gular photon redistr ibution of the line photons is discussed by 
Isenberg et alj dl998bl) . Due to the redistribution of photons in 
angle, the cyclotron line shapes also vary significantly. We have 
fitted the fundamental feature with the continuum spectral func- 
tion multiplied with a Lorentzian in absorption for the absorption 
feature and two Lorentzians in emission for the emission wings. 
Fig.[8]shows characteristic parameters of the fundamental CRSF 
as obtained from these phenomenological fits. 

The line position of the fundamental CRSF varies little with 
the viewing angle and thus cannot account for the amount of 
change in line positions observed for some sources with phase 
or during various observations. Therefore, different explanations 
have to be sought. A possible scenario could be a variation of the 
magnetic field with angle. The line depth increases with cos 9 
in the case of cylinder geometry, and decreases for slab geom- 
etry. This is easily understood from Eqs. (0 and fllOl l, as they 
predict the largest optical depth for small 9 for cylinder, and 
for large 9 for slab geometry. For cylinder geometry, the line 
width increases clearly with cos 9 as expected from the angle- 
dependence of relativistic Doppler broadening (see Eq. [T3] be- 
low) and reinforced by the increasing optical depth with cos 9. 
The decreasing optical depth with angle for slab geometry in- 
stead suppresses the trend of the line broadening with cos 9. 

4.3. Line energies vs. magnetic field strength B 

Line profiles for different magnetic field strengths and other- 
wise fixed parameters are shown in Fig. [9] The figure clearly 
shows the approximately linear progression of the centroid line 
energies towards higher energies with increasing magnetic field 
strength expected from the 12-B-12 rule. The relativistically cor- 



Fig. 9. Line profiles for different values of the magnetic field 
strength. The B-field increases from bottom to top: B/B^t = 
0.05,0.06,0.07,0.08. We assume cylinder geometry, a constant 
temperature kT e = 3keV, an optical depth tj — 3 • 10~ 3 and a 
viewing angle 9 within /j = cos 9 = 0.25. The continuum spec- 
trum has the shape of a power law with photon index a = 2.0 
and with an exponential cutoff at the energy Ef n = 40 keV. 



rect line ratios for a uniform field are implied by Eq. ((3). For 



field are implied by hq Q 
some sources, e.g., V0332+53 dPottschmidt et all [2005h this be- 
havior is confirmed by observational data. As mentioned above, 
there are also sources for which the line energies of the harmon- 
ics are much less harmonically spaced than can be accounted 
for by relativistic effects. One possible explanation is a model 
in which the magnetic field is varying locally within the line- 
forming region. Estimates of the dipole magnetic field of pulsars 
from torque theory and measurements of the surface magnetic 
field from CRSF detections hint at a magnetic field of more 
complex structure than a simple dipole. As the exact nature of 
the magnetic field is unknown, different scenario s for non-dipole 
magnetic field structure have been inves tigated dBlandford et al.L 
1983; lUrpin et al.L Il986t lAronsl 119931) . Besides magnetic field 
gradients which arise due to dipole variations with the height 
of the magnetosphere, surface field va riations can result, e.g ., 



from small-scale crustal field structures (Blandford et al., 



1983 



198 



or thermomagne tic field evolution effects ( Urp in et aU 
iNishimura (120051) recently presented a study of the line ratios for 
a line-forming region of slab geometry threaded by a magnetic 
field which linearly varies with the height. Their approach was 
based on a model bv lGil et all (120021) who assumed the presence 
of a star-centered dipole from a fossil field in the core superposed 
by a crust-anchored dipole anomaly from crustal field struc- 
tures. Applying Feautrier methods to solve the radiative transfer, 
Nishimura (2005) found that the line ratios significantly increase 
if the B-field decreases upwards, and decrease vice versa. We 
have adapted this toy model for a first study of magnetic field 
spread with our Monte Carlo approach. With some simplifying 
assumptions considering the geometry and the angular redistri- 
bution for the case of a non-constant B-field, we can confirm the 
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Fig. 10. Comparison of CRSFs for a uniform and a non-uniform 
magnetic field in a line-forming region of slab 1-0 geometry. 
Bottom: B/B cA = 0.05 -> 0.055, Middle: B/B crit = const. = 
0.055, Top: B/B clit = 0.06 -> 0.055. Results are shown for 
= cos 9 = 0.6. Otherwise the same setting as in Fig. [9] is 
used. The vertical dashed lines mark the line positions obtained 
from a phenomenological fit (continuum multiplied by three 
Lorentzians in absorption for the lines and two Lorentzians in 
emission for the emission wings of the fundamental line) of the 
first three CRSFs. 



t rend o f line-ratio increase and decrease proposed by Nis himural 
(2005). An example is shown in Fig. [TO] where we compare cy- 
clotron lines for a constant, a linearly decreasing, and a linearly 
increasing magnetic field. For the non-constant magnetic field 
we prescribe a linear variation of the field strength in discrete 
steps within the line forming region of up to 10%. The funda- 
mental line appears widely unchanged in shape and position, 
as it is formed in the upper scattering layers, where the non- 
constant B-field was set to have the same value as the constant 
one. Line photons which have been scattered out of the line of 
sight or redistributed in energy in lower layers are replaced by 
spawned photons from scattering in higher layers. By contrast, 
the higher harmonic lines change in position and shape. Here, 
contributions from all layers of different depth are important for 
the final line profile. Absorption features from photons at low 
layers are not refilled. Hence, with the change of the resonant en- 
ergies with the height of the line forming region the lines become 
wider (proportional to the amount of variation in B) and their fi- 
nal centroid energy is shifted. From fitting Lorentzians to the 
first three lines we obtained the line energies and line ratios for 
a constant, increasing and decreasing magnetic field. In Table Q] 
we list the fitted line positions and the fitted line ratios. The the- 
oretical values of the line energies and line ratios after Eq. (0 
are shown for completeness. As expected, these values are sim- 
ilar to the fitted ones for a constant magnetic field, although not 
identical, as the non-Gaussian and non-Lorentzian line shapes 
especially of the fundamental line restrict the fit quality notice- 
ably. 



Table 1. Comparision of line positions and line ratios from the 
fits of the spectra shown in Fig. [TO] The line energies are given in 
keV. For the case of a constant magnetic field, the line energies 
are also calculated from Eq. d3). 
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Fig. 11. Variation of the line shapes with the temperature (kT e = 
5, 10, 15,20keV from bottom to top). Spectra are shown for a 
magnetic field B/B a [ t = 0.07 and otherwise the same parameters 
as in Fig. [9] 



Line ratios of strongly anharmonic nature have been ob- 
served. For 4U 01 15+63, the observed spacing of the line ener- 
gies is smaller than expected, yielding line r atios of (2.8 + 0.05 : 
1.9 + 0.05 : 1) (San tangelo eTall Il999l) or (2.71 + 0.13 : 
1.73 ± 0.08 : 1) dHeindl et all 1 1 99S1) for the first three harmon- 
ics. These deviations are comparable to the ones we obtained for 
an increasing magnetic field in our example. Using a different 
modeling approach, and tun ing his parameters to the case of 4U 
01 15+63, Nishimura (2005) has already shown that these line ra- 
tios can be reproduced by the assumption of a variable magnetic 
field in the line forming region. An example for a source where 
the line ratios may be higher than what is expected f rom Eq. <[3j 
is Vela X-l, where Krevkenbohm et alj dl999ll2002l) have found 
a coupling of the first harmonic energy to the fundamental line 
energy > 2 in RXTE data. 

4.4. Influence of the plasma temperature 

The line width is determined by the energy and angle-dependent 
shape of the scattering cross section and smeared out due to ther- 
mal Doppler broadening. In Fig. Q~T] cyclotron line shapes are 
depicted for varying parallel electron temperature T e , for cylin- 
der geometry, fixed parameters B, cos 6 and for a fixed con- 
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tinuum shape. The hotter the plasma, the wider and the more 
asymmetric become the lines. The width of the lines is due to 
a combination of the natural line width and Doppler broadening 
dHarding & Lail 120061). Doppler broadening gives a Full Width 
Half Maximum of ( Trumper et all Q977; Meszaros & Nagel, 
fl985h 



8 ln(2)fcr e 



cos G\ E c 



(13) 



For increasing cos 9 the line shapes become more asymmetric. 
In our simulation s the plasma temperature is a free param eter. 
From theoretical dLamb et dlll99(tllsenberg et aU [l998b) and 
observational studies dHeindl et all I2004I) . a relation between 
observed magnetic field and temperature has been proposed as 



kT e > 0.2 E c 



(14) 



This would correspond to a temperature for the depicted setting 
(E cyc ~ 30 keV) of at least 6 keV, i.e. somewhere in between the 
bottom spectrum and the second spectrum from the bottom in 
Fig. CD- 



4.5. Continuum shape and photon spawning 

The shape of the cyclotron lines is sensitive to the continuum 
shape. In particular, the fundamental line shape and its emission 
features vary significantly. For better illustration we show the 
case of internally irradiated plasmas where the emission wings 
are strongest. The dependence of the line shapes on the contin- 
uum can be understood when considering the photon redistri- 
bution in energy, especially due to photon spawning. Figure [12] 
shows the change of the line profiles for a flat input continuum 
spectrum, when allowing only for electron transitions between 
the ground Landau state to the first Landau level, or for photon- 
electron scattering leading to up to three harmonics. In the for- 
mer case, a single absorption line forms. The more harmonic 
scatterings are allowed for, the more lines form, while the fun- 
damental and lower harmonics become shallower with growing 
emission wings. Integrating the photon flux only over the en- 
ergy range including just the fundamental line and its emission 
wings (E < 18 • keV B\2), we find that the spawned photons ac- 
count for as much as 34/64/73 % (n < 2, 3,4) of the flux for 
cylinder and for 11/32/43 % of the flux for slab 1-1 geometry. 
Considering the whole energy range, the percentage of spawned 
photons is 25/52/65 % of the total flux for cylinder and 3/15/25 
% of the total flux for slab geometry. Note that these numbers 
are representative of the extreme and fairly unrealistic case of a 
flat input continuum. However, they illustrate well that the line 
shapes should change with the spectral hardness of the incident 
continuum, where harder spectra exhibit more emission features 
near shallower lines. Line profiles for a power law with a photon 
index a — 1 with exponential rolloff at different folding energies 
£foid = 5, 15keV, and for a pure power law continuum spec- 
trum are shown in Fig. [13] In this case, the number of photons 
around the first line increases by a factor of 1.1 and 2.2 for the 
two harder spectra with respect to the softer one. 

5. Applications 

5.1. A new XSPEC model for cyclotron lines: cyclomc 

At present, a physical model for the analysis of cyclotron lines 
is lacking. Typically, for each line independently, multiplica- 
tive phenomenological model components such as Gaussian , 
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Fig. 12. Changes of the line profiles for consideration of one 
to three harmonics for the scattering processes (bottom to top). 
The bottom graph is overplotted in all cases as a dotted line to 
guide the eye. The input continuum spectrum (dashed line) is 
flat and therefore overemphasizes the effects of emission wings. 
Left panel: cylinder geometry; right panel: slab 1-1 geometry. In 
both panels, lcT e = 3keV, cos6» e [0.125,0.250), t t = 3 • 10~ 3 
and B/B CTit = 0.05. 
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Fig. 13. Line profiles for different continua. All continua have 
the general form of a power law with a photon index a = 1 .0 
(top). In the middle and bottom plots, the spectra have an addi- 
tional high energy cutoff at the folding energies £f id = 15keV 
for the middle and Ef n = 5keV for the bottom plot. We 
show the results for cylinder geometry, lcT e = 3 keV, cos e 
[0.375, 0.500), T T = 3 ■ 10~ 3 , and B/B cr « = 0.05. 
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or Lorentzian absorption lines are used. They yield simplified, 
smoothed line shapes; various CRSFs are fitted independently 
from each other one by one. Although it is possible to fit ob- 
served data well with this phenomenological approach, the loss 
of physical information is very unsatisfactory. Overlaying sev- 
eral Gaussians or Lorentzians renders it also difficult to infer the 
magnetic field strength with the line centroid position, and does 
not permit one to distinguish, e.g., thermal and other broadening 
of the CRSFs. The fact that the line positions are independent 
of each other in this approach invokes further information loss. 
With the increasing quality of the energy resolution of todays' 
detectors it is possible to see more details of the line features 
calling for a physical model instead of a phenomenological one. 

For purposes of data analysis, we have implemented a cy- 
clotron line convolution model called cyclomc as a local model 
for fitting CRSFs in XSPEC. Based on the Green's functions of 
Sect. 13.41 cyclomc is designed to fit up to four CRSFs simultane- 
ously. Not only the magnetic field strength and the temperature, 
but also the optical depth and the ratios of line positions, line 
widths and line depths are determined by the underlying physi- 
cal picture. The quality of the fit hence permits conclusions on 
the accuracy of this picture. As shown in Sect. 14.31 the line ra- 
tios could be a sensible indicator for instance for magnetic field 
variations along the line of photon propagation. The model is ob- 
viously restricted to the physical assumptions and the parameter 
input chosen for the Monte Carlo simulations. During fitting, 
it interpolates the results on the assumed input parameter grid 
which has been described in Sect. 13.41 Since cyclomc is a con- 
volution model, self-consistent line shapes are obtained for any 
given continuum. 

5.2. Observability of the line features 

Before applying cyclomc to an example of observational data, 
we want to address the more general question of the observabil- 
ity of the predicted CRSFs in real source data. It is a fundamental 
question to assess which features of the model still appear when 
folded with the detector responses o f todays' observatories . This 
question had already been posed by Isenberg et al. ( 1998b) con- 
sidering the issue of different geometries and densities related 
to the prediction of significant emission wings. However, we 
are not aware of a related study up to date. In Fig. [14] we 
show the simulated spectrum for a source whic h we assign an 
ideal spectrum of the f orm of the npex model (Mihara, 1995; 
Makishi ma et all[T999l) . a negative and positive power law with 
a common high-energy cutoff, 



F(E) = A(E- a > + f ■ E +a2 ) exp 



■— ) 

Efou I 



(15) 



folded with a chosen set of the Green's functions cyclomc for 
slab 1-1 geometry. The flux is appropriate for typical HMXB ob- 
servations. For an assumed observation time of 20 ks this spec- 
trum was folded with the RXTE HEXTE response and back- 
ground was added. We then fitted our fake spectrum, using a 
npex com ponent for the cont inuum and two Gaussian absorp- 
tion lines dCoburn et all 120021) to model the simulated CRSFs. 
Fig. [14] shows this fit as well as the residuals for a fit of the 
continuum component and for a fit of continuum and line com- 
ponents. The residuals from the continuum shape represent our 
simulated line shapes as we have used the same continuum com- 
ponent for the simulated and for the model spectrum. Fitting the 
first two lines with Gaussians, the emission wings stay very pro- 
nounced in the residuals. For our simulated spectrum we expect 




40 50 
Energy [keV] 

Fig. 14. HEXTE simulations. Upper panel: A simulated spec- 
trum of the form npex ■ cyclomc (slab 1-1 geometry, B\2 = 3.0, 
kT e = 5.5 keV, t t = 5 ■ 10~ 4 , cos6» = 0.5, z = 0.3, a x = 1.5, 
a 2 = 1.5, £ fo id = 6.16keV, A = 4.6 ■ 10~ 3 , / = 1.17) is fit- 
ted with a npex multiplied with two Gaussian absorption lines 
(XSPEC: gauabs). Upper panel: data and fitted model. Middle 
panel: residuals for the simulated spectrum fitted with a npex 
continuum. Bottom panel: residuals for the continuum and lines 
fit shown in the upper panel. The strong emission wings of the 
fundamental line in the simulated spectrum are clearly observ- 
able in the residuals. 



to observe the fundamental CRSF at E° bs =26.1 keV (B 12 = 3.0, 
z = 0.3, cos0 = 0.5). The Gaussian fit of the fundamental line 
gives a centroid energy of 24.77 + 0.03 keV instead. This results 
indicates that the asymmetric line shapes could introduce a sys- 
tematic uncertainty in line energy when modeling observed data 
with Gaussian or Lorentzian shapes. We also note that the lines 
seem very prominent even for low values of continuum optical 
depth. 



5.3. Comparison to observational data 

During the recen t outbursts of the t ransie n t sources V0332+53 
and A0535+26 dPottschmidt et all [20051: iKrevkenbohm et all 



120051; iKretschmar et all 120051; iTerada et all 120061) more data 



from accreting X-ray pulsar systems containing strong cyclotron 
resonance scattering features were collected. We have fitted 
RXTE and INTEGRAL data of the January 2005 outburst of 
V0332+53 as a first example for th e comparison of cyclomc 
with real source data. Discovered bv lTanakal (119831) . V0332+53 
was the fo urth accreting X-ray pul sar system where CRSFs 
were found (Makishima et al., 1990a b ). Its third CRSF was dis- 
cover ed in RXTE data dCoburn et all 120051: iPottschmidt et all 
120051) and confirmed by su bsequent INTEGRAL observations 
( Kreykenbohm et al., 2005). In Fig.[15]we show a fit of the first 
two CRSFs in the INTEGRAL data of V0332+53. The fit of the 
RXTE data of V0332+53 during the same outburst is not shown 
as the principal results turned out very similar for both instru- 
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Fig. 15. Fit of the Green's functions CRSF model cyclomc to 
the line features of a phase-averaged spectrum of V0332+53 ob- 
tained from INTEGRAL IBIS observations. The data and fitted 
model (upper panel) along with Ax residuals are shown (three 
lower panels). The residuals are depicted for (a) a continuum 
fit, and for two different continuum and line fits using (b) slab 
1-1 geometry, and (c) slab 1-0 geometry combined with partial 
covering. The continuum shape is taken as a power law with 
a Fermi-Dirac cutoff (fdcut). For the best fit parameters, see 
Tabled 

ments. We modeled the data with a power law with a Fermi- 
Dirac cutoff (fdcut), 



Table 2. Best fit parameters for the data and model shown in 
Fig. Q3] Uncertainties are at the 90% confidence level for one 
interesting parameter; x 2 = 28.4 for 13 d.o.f. 



name [unit] 



Bn [B/(10 L 
kT t [keV] 
r T [10- 3 ] 
cos6> 

z 
a 

£ cut [keV] 
£ foM [keV] 
Ai 
A 2 



G)] 



fit value 

3 05 +u ui 
10 2 +03 

iu -0.3 

3 +0 

0.06 + « 02 

23 +0 02 
-fl!B 
94+o.os 

12 - 8+ p 



7.5 



+0.1 



trum identical to the theoretical one would indeed still exhibit 
detectable emission wings in real data. These strong emission 
features are clearly not present in the spectrum of V0332+53. 
We therefore conclude that the 1-1 geometry of an internally ir- 
radiated plasma is not a completely valid physical assumption. 
Considering the source photons to be injected at the bottom of 
the slab (1-0 geometry) or according to a probably biased distri- 
bution of incident photon production in the plasma seems to be 
more realistic. Second, we see as a general problem with the fit 
that the line depths cyclomc yields for rather low tj are very 
deep compared to the observations. The scenario of partial cov- 
ering of the emergent radiation can solve this issue. In this case, 
only part of the emergent radiation is assumed to pass the region 
of line formation and is reprocessed with cyclomc. A different 
possibility to ensure shallower line depths, however, would be 
the introduction of a magnetic field gradient along the B-field 
vector in a cylindrical line-forming region (O. Nishimura, pri- 
vate communication). Besides yielding shallower lines, such a 
scenario might also improve the fit due to variations in the line 
ratios. 



F(E) = A ■ E~ a 



exp| — 1+1 

£fold 



(16) 



folded with the CRSF Green's functions (cyclomc) for different 
geometries, and assuming 5% systematics for the IBIS data. The 
best fit was obtained for the slab 1-0 geometry in combination 
with the scenario of partial covering and gave a reduced x 2 ed of 
2.2 for 13 degrees of freedom. Fig. [T51 shows the best-fit model 
along with the data and the residuals for fitting the continuum 
component, fitting cyclomc for slab 1-1 geometry, and fitting 
cyclomc for slab 1-0 geometry, combined with partial cover- 
ing, to the data. The best-fit parameters are listed in Table [2] 
While the slab 1-0 geometry is set apart from the other geome- 
tries by exhibiting the weakest emission wings, the scenario of 
partial covering effectively reduces the line depths of the model 
spectrum (compare Sect. 15.21) . We stress that these fits are pre- 
liminary, serving at this stage of the model design as a proof 
of concept only and as an important first step to help us in un- 
derstanding the physics of CRSF formation. The rather high x 2 
should be seen in the light that this is the first time that a realistic 
CRSF model has been fit to real observational data. The fact that 
we can assess the physics via a direct comparison to data enables 
us to draw conclusions concerning the underlying physical pic- 
ture. First of all, the use of the 1-1 geometry results in emission 
wings which are clearly too pronounced for the observed data. 
We have shown in Sect. |5.2| that a simulated source with a spec- 



6. Summary and Discussion 

We have performed an in-depth study of the formation of cy- 
clotron resonance scattering features in the spectra of highly 
magnetized accreting neutron stars. In particular, we have dis- 
cussed the influence of the magnetic field strength, the plasma 
temperature, the angle of radiation and the seed continuum spec- 
tral form onto the shapes of cyclotron line features. 

Our study, being based on the Monte Carlo simulation code 
oflArava & Harding|(ll996lll999l) and lArava-Gochez & Harding! 
(2000), emphasizes several issues which have been pointed out 
by those authors before. Examples are Figs. [4] and |5J where we 
have discussed the optical d epth progression matchi ng our il- 
lustrations to Figs. 4-6 from lArava & Harding! y999) or Fig. Q 
for the study of angular redistribution s i milar to the illustra- 
tions used in Arava-Gochez & Harding (2000). With respect 
to their previous code our model is different in the following 
points: most importantly, we did not restrict ourselves to the 
study of hard continua, but chose a Green's functions approach, 
thus gaining independence from a priori chosen forms of inci- 
dent radiation. We relaxed the geometrical constraints on the 
line-forming region to include the case of a bottom-illuminated 

I I J I 

slab as it has bee n studied by, e.g., Isenb erg et al.l ([1998b) and 
Nishimura (2005). Moreover, the higher resolution in angle and 
energy binning permits to illuminate the complex form of the 
fundamental feature in detail. The time-consuming calculations 



G. Schonherr et al.: A model for cyclotron resonance scattering features 



13 



on a huge parameter grid in principal allow for a systematic com- 
parison of our results to observational data, which is facilitated 
by the implementation of our local model cyclomc for XSPEC. 

Beyond this, we were able to confirm several results from 
other authors obtained with different numerical approaches. In 
Fig. |7] (Sec. 12.11 ) we have shown the angular redistribution of 
photons for cylinder and slab geometry and for two different 
values of the plasma optical depth. The percentage of redis- 
tributed photons per (cosfl)-bin increases with 6 for cylinder 
and decreases with 6 for slab geometry as shown in a sim- 
ilar study of an internally irrad iated line forming region by 
lArava-Gochez & Ha rding (2000). For a higher optical depth and 
for cylinder geometry, how ever, the photon distrib ution flattens, 
similar to results shown bv llsenbergetafl (H998b) for the com- 
parable case of O = n/2 and slab 1-1 geometry. The study of the 
variation of the cyclotron li ne ratios for a no n-uniform magnetic 
field picks up an idea from Nishimura (2005), who investigated 
this case for a fi-field which varies linearly with height, slab ge- 
ometry, and for similar optical depths. We have confirmed the 
trend of an increase of the line ratios with a decreasing magnetic 
field and a decrease of the line ratios with an increasing magnetic 
field within a line-forming region of slab 1-0 geometry (Fig. [TO]). 

In section 12.11 we reported theoretical predictions of the 
model independently from observational data analysis. Except 
for the study of the line ratios, outlined above, all analysis as- 
sumed a uniform magnetic field in the line-forming region. A 
key result is the study of the variation of the line parameters 
of the fundamental feature with the magnetic field (Fig. [9J, op- 
tical depth and angle (Figs. @]-[8j, and temperature (Fig. ITTb . 
which is consistent with predictions from other authors. We do 
not see a significant variation of the fundamental line energy, 
E cyc , with the angle, and thus must exclude this simple sce- 
nario as an explanation for observed phase dependent variations 
of E cyc . Omitting the emission wings, the depth of the funda- 
mental with respect to the continuum flux is rather stable over 
6 whereas the line width varies significantly with 6 for cylinder 
and slab geometry. For both geometries the lines become wider 
towards higher cos 9; for slab geometry an initial decrease for 
cos 6 < 0.25 is observed. The variation of the overall line fea- 
tures for different magnetic field strengths and different temper- 
atures was investigated. Obviously, the positions of the CRSFs 
are directly linked to the B-field strength (see Eqs. Q] and [3]). 
However, Fig. [9] also shows that changes with B as to the line 
shapes are rather insignificant. On the other hand, the line shapes 
vary strongly with increasing temperature, where more asym- 
metric, Doppler-broadened lines arise for higher plasma temper- 
atures. Furthermore, we have studied in depth variations of the 
line shapes with the incident continuum shape. In particular, the 
shape of the fundamental line changes with the continuum shape, 
an effect which can be understood from photon redistribution, 
mainly due to photon spawning in hard continua. As a result, for 
hard spectra for instance the emission wings are much more pro- 
nounced than for softer continua. This dependence of the CRSFs 
on the continuum in principle also allows for conclusions on the 
continuum shape when modeling cyclotron resonance scattering 
features. 

This study aims at meeting the interests of observers in an- 
alyzing cyclotron lines. We therefore chose a broad approach to 
the topic of line formation. Fitting CRSFs in V0332+53 with 
cyclomc indicates that we can indeed assess CRSFs in real ob- 
served source data with a physical model. This is the first time 
such a simultaneous fit of several CRSFs with a realistic, phys- 
ical model has been attempted. As outlined in Sect. 15.31 mod- 
ifications of the underlying physical scenario guided by these 



preliminary fits should help in further improving the compari- 
son of the model and real data. In Sect. 15.21 we have assessed 
the general question of the observability of the theoretically pre- 
dicted line shapes. At the early stages of cyclotron line observa- 
tions it was not clear whether emission features - if they were 
present in the data - wou ld be observed or just sme ared out by 
the det ectors. For instance Isenber g" et alJd!998bl) and[Nishimura 
d2005l) observed that the scenario of a radiation source at the bot- 
tom of a slab as geometry for the line forming region leads to 
less emission features in the spectra than an internally irradiated 
plasma. We have shown in section I3T21 that such strong emission 
features as predicted by our scenario should indeed be observ- 
able by the instruments on todays' observatories. The fit results 
from V0332+53 data demonstrate that large emission wings can 
be ruled out for this spectrum. 

This work is ongoing, aiming at further generalizations of 
the CRSF model. We int end to permit angular anis o tropy of the 
continuum photon flux. Araya-Gochez & Harding (2000) have 
already studied this case for special angular distributions like 
peaked emission. We strive for the calculation of our Green's 
functions independently not only of the continuum energy but 
also of the continuum angular distribution to generalize our 
model to arbitrary angular distributions of the incident pho- 
tons. The first issue to be investigated further, however, will be 
the observed discrepancy in theoretical and observational line 
depths. More realistic B-field gradients could account for shal- 
lower lines. We have shown that the consideration of a scenario 
of partial covering, which effectively reduces the line depths, 
significantly improves the fit quality. We therefore think that 
the consideration of non-constant magnetic fields within the line 
forming region will play a major role in the process of better 
modeling and understanding cyclotron lines. This issue will be 
subject of a forthcoming paper once having realized a systematic 
comparison of cyclomc for different scenarios for line forma- 
tion to a larger set of observational data. 
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